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We present a scheme to efficiently simulate, with a classical computer, the dynamics of multipartite 
quantum systems on which the amount of entanglement (or of correlations in the case of mixed-state 
dynamics) is conveniently restricted. The evolution of a pure state of n qubits can be simulated by 
using computational resources that grow linearly in n and exponentially in the entanglement. We 
show that a pure-state quantum computation can only yield an exponential speed-up with respect 
to classical computations if the entanglement increases with the size n of the computation, and gives 
a lower bound on the required growth. 
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In quantum computation, the evolution of a multipar- 
tite quantum system is used to efficiently perform com- 
putational tasks that are believed to be intractable with 
a classical computer. For instance, provided a series 
of severe technological difficulties are overcome, Shor's 
quantum algorithm Q can be used to decompose a large 
number into its prime factors efficiently — that is, expo- 
nentially faster than with any known classical algorithm. 

While it is not yet clear what physical resources are 
responsible for such suspected quantum computational 
speed-ups, a central observation, as discussed by Feyn- 
man |2( , is that simulating quantum systems by classical 
means appears to be hard. Suppose we want to simulate 
the joint evolution of n interacting spin systems, each 
one described by a two-dimensional Hilbcrt space H.2- 
Expressing the most general pure state G TL 2 ® n of 
the n spins already requires specifying about 2™ complex 
numbers c^...^ , 
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where {|0), |1) € TL2} denotes a single-spin orthonormal 
basis; and computing its evolution in time is not any 
simpler. This exponential overhead of classical computa- 
tional resources -as compared to the quantum resources 
needed to directly implement the physical evolution by 
using n spin systems- strongly suggests that quantum 
systems are indeed computationally more powerful than 
classical ones. 

On the other hand, some specific quantum evolutions 
can be efficiently simulated by a classical computer - 
and therefore cannot yield an exponential computational 
speed-up. Examples include a system of fermions with 
only quadratic interactions [j| , or a set of two- level sys- 
tems or qubits initially prepared in a computational-basis 
state and acted upon by gates from the Clifford group Q . 
Recently, Jozsa and Linden || have also shown how to 
efficiently simulate any quantum evolution of an n-qubit 
system when its state factors, at all times, into a prod- 
uct of states each one involving, at most, a constant (i.e. 
independent of n) number of qubits. 



Here we show how to efficiently simulate, with a classi- 
cal computer, pure-state quantum dynamics of n entan- 
gled qubits, whenever only a restricted amount of entan- 
glement is present in the system. It follows that entan- 
glement is a necessary resource in (pure-state) quantum 
computational speed-ups. More generally, we establish 
an upper bound, in terms of the amount of entangle- 
ment, for the maximal speed-up a quantum computation 
can achieve. An analogous upper bound, but in terms of 
correlations (either classical or quantum), also applies to 
quantum computations with mixed states. 

For simplicity sake the analysis is focused on a compu- 
tation in the quantum circuit model. Thus we consider 
a discretized evolution of the n qubits, initially in state 
10)®", according to a sequence of poly(n) (i.e., a num- 
ber polynomial in n) single-qubit and two-qubit gates. 
We recall, however, that any evolution of n qubits ac- 
cording to single-qubit and two-qubit Hamiltonians can 
be efficiently approximated, with arbitrary accuracy, by 
the above circuit model, so that the present results also 
apply to this more general setting [f| . 

Consider, as in Eq. a pure state |^) G H^ n of an 
?i-qubit system. Let A denote a subset of the n qubits 
and B the rest of them. The Schmidt decomposition SD 
of I \P) with respect to the partition A:B reads 

\9) = £\ a \*£ ] )®\* l ? i ), (2) 

a=l 

where the vector \<f> [ a ] ) (|$L B1 )) is an eigenvector with 
eigenvalue |A a | 2 > of the reduced density matrix p 1 ^ 
(p[ s l), whereas the coefficient A Q follows from the rela- 
tion ($L A1 |*) = A a |$L S1 ). The Schmidt rank X A is a 
natural measure of the entanglement between the qubits 
in A and those in B 0. Accordingly, we quantify the 
entanglement of state \9) by X, 

X = maxXA, (3) 

A 

that is, by the maximal Schmidt rank over all possible bi- 
partite splittings A:B of the n qubits. We shall say that 
|\&) is only slightly entangled if X is "small". In particu- 
lar, here we are interested in sequences of states {|'3/ T1 )]- of 
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an increasing number n of qubits (corresponding, say, to 
quantum computations with increasingly large inputs). 
In such a context we consider X to be "small" if it grows 
at most polynomially with n, X n = poly(n) ||. 

Definition.— A pure-state quantum evolution is 
slightly entangled if, at all times t, the state ^(i)) of the 
system is slightly entangled — that is, if x(t) is small. A 
sequence of evolutions with an increasingly large number 
n of qubits is slightly entangled if X n (t) is upper bounded 
by poly(n). 

The key ingredient of our simulation protocol is a local 
decomposition of the state |^) 6 Hz® 71 in terms of n 
tensors {rW}f =1 and n-1 vectors {A^}^, denoted 

l#) < y r [1 U [ % [2 U [2] • • • r [ ' ] • • • A [n_11 rH (4) 

Here, tensor is assigned to qubit I and has (at most) 
three indices, T^, , where a, a' = 1, ■ • ■ , X and i = 0, 1, 
whereas A^ is a vector whose components A^l store the 
Schmidt coefficients of the splitting [1 •••/]:[(/ + 1) •• -n]. 
More explicitly, we have 9] 

r \" p[%\[(]r[2]<2 A[2]...rH>» (k\ 

Ql ,- ■ ■ ,a n —i 

so that the 2™ coefficients c^...^ are expressed in terms of 
about (2x 2 + X)n parameters, a number that grows only 
linearly in n for a fixed value of X- This decomposition 
is local in that, as we shall see, when a two-qubit gate is 
applied to qubits I and l+l, only T^, A^ and r^ +1 l need 
be updated. 

Decomposition (@J) (but not X) depends on the par- 
ticular way qubits have been ordered from 1 to n, and 
essentially consists of a concatenation of n—l SDs. We 
first compute the SD of \^f) according to the bipartite 
splitting of the systems into qubit 1 and the n—l re- 
maining qubits [from now on we omit the tensor product 
symbol] , 

I*) = £ A «>I $ £"" ] > (6) 

Ql 

= Er^^lKi)!^ 2 ;-" 1 ), (7) 

i 1 ,a 1 

where in the last line we have expanded each Schmidt 
vector |$a|) = Yli To]* 1 \h) in terms of the basis vectors 
{|0), |1)} for qubit 1. We then proceed according to the 
following three steps: (i) first we expand each Schmidt 
vector |<&q in a local basis for qubit 2, 

l^HENl^'r 1 ); (8) 

(ii) then we write each (possibly unnormalized) vec- 
tor It^'™^) in terms of the at most X Schmidt vectors 
{|$a 2 n ^)}a 2 =i (i- e -; the eigenvectors of pi 3 '""!) and the 
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corresponding Schmidt coefficients X a2 , 

i^: , } = E r S 2 A Si*'" 1 ); ( 9 ) 

a 2 



(Hi) finally we substitute Eq. © in Eq. (|SJ) and the 
latter in Eq. J7J to obtain 

l*>= E rm <i AL 1 lrP^AM| il i 3 )|$L 3 r B] >. (10) 

Iterating steps (i)-(iii) for the Schmidt vectors 

|$L 3 2 " nI ), |$L 4 3 " nI ), • ■ • , \^SLx), one can express state |*) 
in terms of tensors and A^, as in Eq. Q. 

A useful feature of description Q is that it readily 
gives the SD of |W) according to the bipartite splitting 
[l-..Z]:[(I+l)...n], 

i*) = E A Si^r i] )i^r- n] ). (ii) 

eti 

Indeed, it can be checked by induction over I that 

|$£r fl >«— rMAM.-.A^rM, (12) 
meaning that 

= E r *\ hx *l ■ ■ ■ r S!U (i3) 

ai,- -, a i-i 

whereas by construction we already had that 

| $ K|+i)-»] ) r R+i] A [/+i] . . . A [n-i] r N (14) 

which stands for 
1$^ )-»]) = 53 rlWfr • • • Afc^N- | iw . . . in) .(i5) 

The following lemmas explain how to update the de- 
scription of state \^>) when a single-qubit gate or a two- 
qubit gate (acting on consecutive qubits) is applied to 
the system. Remarkably, the computational cost of the 
updating is independent of the number n of qubits, and 
only grows in X as a polynomial of low degree. 

Lemma 1.— Updating the description of state |\I>) 
after a unitary operation U acts on qubit I does only 
involve transforming . The incurred computational 
cost is of 0(X 2 ) basic operations. 

Proof.— In the SD according to the splitting [1 • • • (I— 
1)] : [/■••7i], a unitary operation U on qubit I does 
not modify the Schmidt vectors for part [1 • • • (I — 1)] 
and therefore and A^ (1 < j < I - 1) remain the 
same. Similarly, by considering the SD for the splitting 
[1 • • -I] : [(l+l) • • ■ n], we conclude that also 1*1 and A^ 1 ! 
(' + 1 < j < n) remain unaffected. Instead, changes 
according to 

r'g=E^FS' Va,/3 = 1,.-.,X. (16) 

3=0,1 

Lemma 2. Updating the description |Q of state |\I>) 
after a unitary operation V acts on qubits I and I + 1 
does only involve transforming T^, A^ and T^ l+1 \ This 
can be achieved with 0(X 3 ) basic operations. 
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Proof.— In order to ease the notation we regard \$f) 
as belonging to only 4 subsystems, 



H = J ®Uc®Ud®K>- 



(17) 



Here, J is spanned by the X eigenvectors of the reduced 
density matrix 

^-(«)] = £Ka|, \a) ^ Xl^t il - 1)] h (18) 



and, similarly, K, is spanned by the X eigenvectors of the 
reduced density matrix 



AW)- 



= 5>XtI, |7)-A^|$r 2 >'-]>; (19) 



whereas He and Hd correspond, respectively, to qubits 
I and l + l. In this notation we have 

l*>= E ErS%r^|mi 7 ), (20) 

q,/3,7— 1 

and, reasoning as in the proof of lemma 1, when apply- 
ing unitary V to qubits C and D we need only update 
r[°], A.T^l. We can expand |*') = V|#) as 



where 



X i 

a, 7=1 i,j=0 



U t> 7 — / j / j v kl 1 uB A P l 0f ■ 



B kl 

By diagonalizing p'^ DK \ 



P 



/[DfC] 



(21) 



(22) 



(23) 



= E E< a i a > e ^( X<)* LnXjVI, 

jj'nn' \ ) 

we obtain its eigenvectors {l^'g ^}} 5 which we can ex- 
pand in terms of { | JT) } t° obtain r'l 15 !, 

i*T ] )=E r ri^)- ( 24 ) 

The eigenvectors of p'^ c ^ and A' follow then from 

= E ( r 57 b ')* ^<7l7)M, (26) 



and by expanding each | $ i ) , 



(27) 



we also obtain f^. All the above manipulations can 
be performed by storing 0(x 2 ) coefficients and require 
0(X 3 ) basic operations. 



We now state our main results. We consider a pure- 
state quantum computation using n qubits, and consist- 
ing of poly(n) one- and two-qubit gates and a final lo- 
cal measurement. The simulation protocol works as fol- 
lows. We use tensors r'" and A^ to store the initial 
state \0)® n and update its description as the gates are 
applied ^(J- Recall that in description (@J each qubit 
has been associated a position from 1 to n. In order to 
update \^>) according to a two-qubit gate between non- 
consecutive qubits C and D, we will first simulate 0(n) 
swap gates between adjacent qubits to bring C and D to- 
gether. Computing the expectation value for any product 
operator (e.g. a projection corresponding to a local mea- 
surement) from {r^, A^} is straightforward and can also 
be done with n poly(x) operations. 

Theorem 1. — If through a pure-state quantum com- 
putation X n is upper bounded by poly(n), then the com- 
putation can be classically simulated with poly(n) mem- 
ory space and computational time. 

Theorem 2. — If X n grows subexponentially in n, 
then the quantum computation can be classically sim- 
ulated with subexp(n) memory space and computational 
time. 

Thus, theorem 1 provides us with a sufficient condi- 
tion for the efficient classical simulation of a quantum 
computation, which by extension also applies to generic 
pure-state, multi-particle unitary dynamics generated by 
local interactions |6j . In turn theorem 2 provides us with 
a more general condition under which a quantum compu- 
tation cannot yield an exponential speed-up with respect 
to classical computations. Both theorems follow straight- 
forwardly from the previous lemmas and considerations. 

The above results establish a clear connection between 
the amount of entanglement in a multipartite system and 
the computational cost of simulating the system with a 
classical computer. This suggests a new approach to the 
study of multipartite entanglement, based on the com- 
plexity of describing and simulating quantum systems. 
We propose to quantify the entanglement of a pure state 
| "J/) through measures that indicate how difficult it is to 
express \ty) in terms of local states or, relatedly, to ac- 
count for a local change in the system. An example of 
such entanglement measures is the function 



log 2 X, 



(28) 



which, apart from serving the purposes, has a series 
of other appealing properties: (i) Ex only vanishes for 
product (i.e., unentangled) vectors; (ii) E\ is additive 
under tensor products, E x (^ *') = E x (^>) + E x {^')] 
(Hi) E X monotonically decreases under (both determin- 
istic and stochastic) LOCC manipulations of the system. 
We also note that Ex(*&) is not a continuous function of 
|\&) with respect any reasonable distance [ll| . 

We can rephrase the results of this paper in terms of 
Ex- Notice that the maximum value of in a sys- 
tem of n particles is linear in n. Theorem 1 states that 
an efficient simulation of quantum dynamics is possible 
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whenever Ey/ grows at most logarithmically in n. More 
generally, we have shown how a state \ty) can be given a 
description in terms of local states by using a number of 
parameters that grows linearly in the number of systems 
and exponentially in the amount of entanglement Ex, 

local description ^ nexp(i?x) / 00 n 
of an n-qubit state ~ parameters. 

This expression implies an upper bound, in terms of the 
entanglement, for the computational speed-up a quantum 
evolution can achieve with respect to classical computa- 
tions. 

So far we have only considered pure-state dynamics. 
But if the n qubits are in a mixed state p £ B{Tl2® n ), 
we can regard density matrices as vectors in the space 
of linear operators. By using product expansions and 
the Schmidt decomposition in this space, one can readily 
re-derive the above results, but with the former role of 
entanglement played now by both quantum and classi- 
cal correlations. Thus, an efficient simulation is possi- 
ble if the total amount of correlations (as measured by 
the analog of X) is sufficiently restricted. In particular, 
this results do not rule out the possibility of obtaining 
a computational speed-up through a quantum computa- 
tion with very noisy mixed states 

Finally, [a simple modification of] the simulation proto- 
col discussed in this paper may find practical applications 
as a tool to study quantum systems 0. The results 
of suggest that, at zero temperature, non-critical 
spin-chains typically meet sufficient conditions for an effi- 
cient classical simulation. Perhaps, then, understanding 
the structure of multipartite entanglement is the key to 
achieve efficient simulation of certain multipartite quan- 
tum phenomena. 
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